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Abstract 

We consider SU(A^) QCD^^^ coupled to massless adjoint Majorana fermions, 
where N is finite but arbitrary. We examine the spectrum for various values of N, 
paying particular attention to the formation of multi-particle states, which were 
recently identified by Gross, Hashimoto and Klebanov in the N = oo limit of 
the theory. It is believed that in the limit of vanishing fermion mass, there is a 
transition from confinement to screening in which string-like states made out of 
adjoint fermion bits dissociate into stable constituent "single particles". In this 
work, we provide numerical evidence that such a transition into stable constituent 
particles occurs not only at large N, but for any finite value of N. In addition, 
we discuss certain issues concerning the "topological" properties exhibited by the 
DLCQ spectrum. 
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1 Introduction 



Solving for the non-perturbative properties of physically realistic quantum gauge theories 
is typically an intractable problem. In order to gain some insights, however, a number of 
lower dimensional models have been investigated in the large (or planar) approxima- 
tion, with a plethora of examples emerging in the last few years (see reference for an 
extensive review). 

In this work, we will consider the 1 + 1 dimensional SU(A^) gauge theory of QCD 
coupled to massless adjoint Majorana fermions, which is believed to exhibit the property 
of screening 0, [T^. A similar theory with complex adjoint fermions has also been con- 
sidered [|TT], [T^. It will be advantageous to quantize the theory on the light-cone, and to 
adopt the light-cone gauge. It is then a straightforward task to extract numerical bound 
state solutions via an application of Discrete Light-Cone Quantization (DLCQ) 

Various non-perturbative studies of this model already exist in the literature [|, ^, 0, 
but in more recent work by Gross, Hashimoto and Klebanov (hereafter 'GHK'), it was 
suggested that in the massless fermion limit, the spectrum becomes continuous above a 
certain threshold. This was supported by the presence of states in the spectrum that have 
a mass consistent with the dynamics of two freely interacting stable particles. In this 
context, identifying either "single particle" or "multi-particle" states requires a careful 
analysis of the spectrum, since an analysis of the explicit Fock state content does not 
distinguish these states unambiguously. 

It is worth clarifying this last remark to avoid possible confusion. Firstly, since the 
limit = cxD is assumed in GHK, Fock states are single traces of fermion creation 
operators, and correspond intuitively to single closed strings of adjoint fermion bits. 
Multi-trace states (corresponding to multi-string states) do not appear in the analysis, 
since interactions with them are suppressed (the factor 1/A^ plays the role of a string 
coupling constant). Naively, one would view single-trace states as "single particles" and 
multi-trace states as "multi-particles". Although the latter is expected to follow from 
the usual 1/A^ counting, GHK have shown that there may be exceptions to the former; 
namely, states that are superpositions of single-trace Fock states may exhibit masses that 
are consistent with the dynamics of two freely interacting stable particles. 

This property of the large A^ spectrum was perhaps anticipated by the work of Kutasov 
and Schwimmer namely, the presence of multi-particle states in the N = oo 

theory becomes manifest after an appropriate rearrangement of the Hilbert space into 
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Kac Moody current blocks. At finite A^, however, interactions prevent us from 
identifying obvious multi-particle candidates in the Hilbert space[], and we can no longer 
appeal to the correspondence proposed in to argue for the existence of multi-particle 
states. Nevertheless, the formation of multi-particles in the finite spectrum may 
be inferred from the screening properties of the theory [^], since the theory screens for 
any finite A^[@]. It is therefore of interest to determine whether multi-particles actually 
appear in the finite A^ spectrum or not, and we will devote ourselves towards answering 
this question via a numerical study of the DLCQ bound state equations. We also attempt 
to shed light on certain "topological" features exhibited by the DLCQ spectrum. 

Of course, working at finite A^ has a price; the Fock space now admits multi-trace 
states, and the complexity of the bound state problem is dramatically increased. Nev- 
ertheless, we find that the numerical bound state problem is still tractable provided the 
discretization of the light-cone momentum is not too fine, but fine enough to resolve 
certain features of the theory we are interested in. 

The organization of the paper may be summarized as follows; in Section ^ we briefiy 
review the 1 + 1 SU(A^) gauge theory coupled to adjoint Majorana fermions, which we 
formulate in light-cone coordinates. We also discuss features of the DLCQ formulation 
at finite A^ that differ from the large A^ formulation originally discussed in Q . In Section 
^ we tabulate the results of our finite A^ numerical analysis, and compare spectra of 
candidate multi-particle states with mass predictions given by free-body kinematics. We 
conclude our investigation with a summary and discussion in Section ^. 

2 Light-Cone Quantization and DLCQ at Finite 

The action for QCD^^^ coupled to a single Majorana fermion transforming in the adjoint 
representation of U(A^) or SU(A^) is given byQ 



S = d xTy 



1 

47 



:ii 



where the covariant derivative is defined by D^"^ = 9^\E' + i[A^, and the fermion field 

X 



v]/ _ 2 j ^ ) is a two component spinor, each component representing an A^ x A^ 



^For large N, multi-trace states are obvious candidates for multi-particles. 

^We consider V{N) as well as SU(iV), since it was shown in earlier work that the SU(A^) spectrum 
may be obtained by solving for the U(A^) spectrum, and then 'factoring out' U(l) states. Numerically, 
this method is considerably more efficient than solving for the SU(A^) spectrum directly. 
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Hermitian matrix of Grassmann variables. These matrices are traceless if the gauge group 
is SU(A^). The hght-cone quantization of this theory in the hght-cone gauge = was 
carried out in and we refer the reader to that source for details and notation. In the 
present context, we simply note that the light-cone momentum P'^ and Hamiltonian 
may be expressed in terms of Fourier oscillator modes as follows: 



P^ 
P 



dk k blj{k)bij{k), 



dk 



(2) 
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For the gauge group U{N), the fermionic oscillator modes above satisfy the following 
anti-commutation relations 

{k,{k^),bl{k'-)} = S{k^ -~k^)SuS,k. (8) 
If the gauge group is SU(A^) we need to adopt the following set of relations - 



{%(A:+),4(fc+)} = 6ik^ - k^){5u5,u - -^5,^ 



(9) 



and discard the j^b\^{k)b.jj{k) term in P~, since the SU(A^) fermion fields are traceless. 

A number of comments are in order. If we compare the light-cone Hamiltonian P~ of 
equation - which is valid for any - to the large N expression given in |Q , we notice 
that there is an additional bb bb term in with amplitude B{ki). The color index 
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structure of this termQ implies that it is suppressed by since it sphts a single-trace 
Fock state into multi-trace states. Nevertheless, the Z2 symmetry T : hij —* bji of the 
large N theory ^ is also manifest in the finite N formulation, since this additional term 
is easily shown to be invariant under T. 

The D term in equation (^) is obtained by normal ordering quartic terms in P~ , 
and diverges linearly. However, this divergence is entirely absorbed by the Coulomb 
divergence generated by the A term. So the theory is manifestly finite [Q. We also note 
that the D term used in the large N expression for P~ is unchanged if is finite - the 
contributions obtained from normal ordering quartic terms in accordance with the 
SU(A^) relations (||) simply sum to zero. 

In order to implement the DLCQ formulation of the bound state problem [0, we 
simply restrict the light-cone momentum variables ki appearing in equations (0),(0) for 

to the following set of discretized momenta: {-7^, ^j^, ^j^, ■ ■ ■}', i.e. only ode? positive 
integer multiples of P~^/K are allowed, which is equivalent to imposing anti-periodic 
boundary conditions^ for the fermion fields: ipij{x~) = —ipij{x~ + 2ttR). The integer K 
is called the harmonic resolution, and 1/K measures the coarseness of our discretization. 
Physically, 1/K represents the smallest unit of longitudinal momentum fraction allowed 
for each parton. As soon as we implement the DLCQ procedure, which is specified 
unambiguously by the harmonic resolution K, the integrals appearing in the definitions 
for are replaced by finite sums, and the eigen-equation 2P~^P~\'$) = M^|\l/) 
is reduced to a finite matrix problem. Continuum values are obtained by extrapolating 
results to the K = 00 limit. Typically, a computer program is used to generate and 
diagonalize the DLCQ matrix to solve for the mass eigenvalues M^. In the present 
work, we are able to perform numerical diagonalizations for values of K in the range 
3 < K < 16 with the help of Mathematica and a desktop PC. At K = 16, the DLCQ 
matrix has dimensions 375 x 375. A similar finite A^ analysis of a two dimensional 
super symmetric matrix model was performed recently by the authors |p. 



■^Repeated indices are always summed from 1 to iV. 
Periodic boundary conditions are also possible, but convergence in numerical calculations is slower 

i- 
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3 DLCQ Bound State Solutions 



In this section we present the resuhs of our numerical diagonahzations of the DLCQ 
matrix for = 2P+p- for values of K in the range 3 < fsT < 16, and foif| N = 3, 10, 100 
and 1000. Strictly speaking, there is no upper limit on the size of A^, since it appears 
as an algebraic variable in the DLCQ matrix. The mass m of the fermion (see equation 
(0)) is set to zero. The results we obtain for = 1000 agree with the large N results 
presented in ^ to at least six significant figures. 

As we stated earlier, our main objective is to determine whether the multi-particle 
states of the N = oo spectrum also persist at any finite A^. Towards this end, we 
first calculate the masses of various "single particle" states - namely, the two lightest 
fermion and boson single-particle states - for various values of K and A^. The results^ 
are presented in Table ^ and Table 0. 

Note that for small values of K there is no dependence on A^ (after expressing 
in units q'^N/tc). Moreover, for larger values of K, the A^ dependence is surprisingly 
small. This seems puzzling at first, since we know interactions between single and multi- 
trace states are governed by the 'string coupling' 1/A^, and so for A^ = 3, one expects 
considerable multi-trace contributions in a state such as \F2), which turns out to be a 
superposition of mainly five-parton single trace Fock states. Evidently a rather remark- 
able cancellation seems to be responsible for keeping the masses relatively independent 
of A^. In this case, the limit N = oo provides an excellent approximation to the A^ = 3 
case. 

A crucial observation made by GHK |^ in their recent work suggested that many 
of the remaining states in the spectrum are adequately described as two or more freely 
interacting single particles: \F1) ® \F1), \F1) (g) \F2), \F1) l-Bl), and so on. This is why 
the states in Tables |I| and ^ are referred to as "single particles" , since they themselves 
are not seen to decompose into more fundamental stable particles. 

The validity of this scheme may be tested in the DLCQ framework as follows For 
a given DLCQ resolution K, the mass Mpif^Fi{K) of the composite state \F1) (E) \F1) is 

^ We exclude N = 2, since in this case one has to calculate the norm of each Fock state explicitly to 
avoid overlapping states, which is computationally very intensive. 

^ We have also identified a number of other single particle states but they are not relevant for the 
discussion presented here. 
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Fermion "Single Particle" Masses (in units g'^N/ir) 


|F1) 


|F2) 


K 


= 3 


iV = 10 


N = 100 


N = 1000 


K 


= 3 


iV = 10 


N = 100 


N = 1000 


3 


4.5 


4.5 


4.5 


4.5 


3 










5 


5 


5 


5 


5 


5 


12.5 


12.5 


12.5 


12.5 


7 


5.22272 


5.22272 


5.22272 


5.22272 


7 


14 


14 


14 


14 


9 


5.34559 


5.34559 


5.34559 


5.34559 


9 


14.7645 


14.7645 


14.7645 


14.7645 


11 


5.42226 


5.42224 


5.42224 


5.42224 


11 


15.2575 


15.2575 


15.2575 


15.2575 


13 


5.47410 


5.47406 


5.47406 


5.47406 


13 


15.5909 


15.5908 


15.5908 


15.5908 


15 


5.51121 


5.51115 


5.51115 


5.51114 


15 


15.8314 


15.8311 


15.8311 


15.8311 



Table 1: The masses (in units g'^N/n) of the two hghtest "single particle" fermion 
states, and \F2) respectively. 



Boson "Single Particle" Masses (in units g'^N/i:) 


|B1) 


|B2) 


K 


N = 2, 


N = IQ 


N = 100 


N = 1000 


K 


N = 2, 


iV = 10 


N = 100 


N = 1000 


4 


8 


8 


8 


8 


4 










6 


9 


9 


9 


9 


9 


18 


18 


18 


18 


8 


9.49097 


9.49097 


9.49097 


9.49097 


8 


20 


20 


20 


20 


10 


9.78145 


9.78145 


9.78145 


9.78145 


10 


21.2117 


21.2117 


21.2117 


21.2117 


12 


9.97108 


9.97103 


9.97102 


9.97102 


12 


22.0078 


22.0078 


22.0078 


22.0078 


14 


10.1035 


10.1034 


10.1034 


10.1034 


14 


22.5681 


22.5680 


22.5680 


22.5680 


16 


10.2006 


10.2004 


10.2004 


10.2004 


16 


22.9812 


22.9811 


22.9811 


22.9811 



Table 2: The masses (in units q'^N/tt) of the two lightest "single particle" boson 
states, l-Bl) and \B2) respectively. 



given by free body kinematics according to the relation 

Mli^^i(ir) _ M|i(n) ^ MUK-n) ^ ^^^^ 
K n K — n ^ 

where n is any positive integer less than and where Mpiin) and Mfi{K — n) are 
masses of the \F1) particle carrying n and K — n momentum units respectively. In Table 
H, we have applied equation (|T0|) to determine the mass of \F1) ®\F1) for various choices 
of n and K. Since the mass of \F1) for each K is essentially independent of at five 
significant figures, the masses given in Table ^ are applicable for any A^. In Table ^ we list 
the actual masses for bosons that are observed in the DLCQ spectrum for various choices 
of K and A^. Comparing Tables |^ and ^, we see for A^ = 1000 (i.e. essentially 'large 
A^'), that the masses corresponding to two free \F1) particles appear in the actual DLCQ 
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Calculated Masses for (g) 


K 




8 


20 






10 


(20.00) 


22.461 




12 


20.953 


25.127 




14 


(20.891) 


22.315 


27.901 


16 


21.441 


23.887 


30.737 



Table 3: Predicted values for (in units q'^N/tt) for the composite particle \F1) (g) \F1) 
according to the two-free-body formula (p!OD, for various choices of n and K. The different 
numbers appearing in a given row correspond to varying the value of n in equation 
(p!0|). Numbers that are in parentheses correspond to a pair of identical \F1) particles 
(i.e. carrying the same momentum), and are therefore expected to be absent from the 
spectrum because of Fermi statistics. 



Observed Boson Masses M"^ in DLCQ Spectrum 


K 


N 


(in units g'^N/ir) 


8 


3 


18.7312 


20(2) 










8 


1000 


18.7312 


20(2) 










10 


3 






21.1213 


22.461(2) 






10 


1000 






21.1213 


22.461(2) 






12 


3 


20.5603 


20.9917^2) 


23.954 


25.0849(2) 






12 


1000 


20.5048 


20.9532^2) 


24.0033 


25.1275(2) 






14 


3 






21.7082 


22.3756(2) 


26.8465 


27.8003(2) 


14 


1000 






21.6346 


22.3154(2) 


26.9517 


27.9010(2) 


16 


3 


21.2700 


21.4620^2) 


23.2368 


23.9392(2) 


29.7486 


30.5640(2) 


16 


1000 


21.2303 


21.4409^2) 


23.1910 


23.8869(2) 


29.9211 


30.7373(2) 



Table 4: Actual values for the mass squared M'^ (in units g'^N/ir) of several bosonic states 
for different N and K (excluding single particle states). Numbers with a superscript ('^•' 
correspond to an exact n-fold degeneracy in the spectrum. 



spectrum with an exact two-fold degeneracy. In addition, the masses in Table ^ that 
correspond to two identical \F1) particles (in parentheses) are absent from the DLCQ 
spectrum, which is of course consistent from the Fermi statistics of two freely interacting 
identical fermions. The two-fold degeneracy we see here in the bosonic spectrum does not 
occur in the N = oo analysis of GHK , and this can be easily understood by analyzing 
the Fock state content of each doublet. What we find for any doublet at = 1000 is 
that one state is a superposition of essentially single-trace Fock states, while the other is 
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mainly a superposition of two-trace Fock states. Since the Hilbert space for the N = oo 
theory is generated from single-trace Fock states only, the multi-trace states we see here 
will be absent in the N = oo spectrum. 

One remarkable property of the theory, which can be observed from Table |^, is that 
for small values of - say, = 3 (where interactions can no longer be neglected) - 
the exact two-fold degeneracy still survives, and the actual masses deviate only slightly 
(if at all) from the large N results. Nevertheless, the Fock state content of each state 
in a doublet is drastically altered if we vary A^; single-trace and multi-trace Fock states 
now contribute equally in each state for small A^. Evidently, significant cancellations 
between 1/A^ contributions in the finite A^ Hamiltonian must be occurring in order to 
keep masses relatively independent of A^. We should point out that the additional term 
in equation (^) for P~ that distinguishes it from the N = oo Hamiltonian is crucial in 
these calculations - no such mass degeneracy would be observed if it was omitted. It is 
still unclear at this point whether the very small deviations from the predicted masses 
in Table |^ and the mass of the doublets appearing in the finite A^ DLCQ spectrum in 
Table ^ will remain in the continuum K ^ oo limit. There is already evidence in Table ^ 
suggesting that the deviations are diminishing for larger values of but clarifying this 
issue will certainly require solving the DLCQ spectrum for larger values of and we 
leave this for future work. 

As was remarked in the work of GHK p|, the remaining boson states listed in Table 
^ agree with the masses presented in Table ^ up to 1/ K] i.e. as we increase the harmonic 
resolution, the agreement becomes more exact. These states are therefore multi-particle 
states that give rise to two-body continua in the continuum limit K —>■ oo. We also note 
that this picture is not disturbed for relatively small values of A^ (such as A^ = 3), and 
so the observations made by GHK in the context of the N = oo theory also appear to be 
valid at any finite A^. In particular, for any finite A^, our results are consistent with the 
interpretation that the spectrum of the (continuum) theory becomes continuous above 
the threshold mass = 4M|,i, where Mpi is the mass of the lightest stable particle in 
the continuum theory. 

An identical analysis may be performed for the composite state \F1) ® \F2), and we 
find once again that the masses predicted by equation (^) emerge in the actual DLCQ 
spectrum as mass doublets for any finite A^. Agreement is more precise for large values 
of A^, but we nevertheless obtain good accuracy even for A^ = 3 as we did for \F1) \ F1) 
in Table ^. For composite particles with fermion statistics, such as \F1) (g) \B1), the 
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masses predicted by equation (p!0| ) also emerge in the DLCQ spectrum for any A^, but 
we no longer see the exact two-fold degeneracy that was observed in the bosonic sector. 
The fermion states that exactly satisfy equation (^) for large are not observed in the 
N = oo analysis of GHK, since these states are a superposition of two-trace Fock states 
(one fermion, and one boson) for large values of A^. For = 3, such states involve a 
superposition of both single and multi-trace states. 

As was pointed out in a large portion of the "multi-particle" spectrum in the 
fermionic sector can be summarized surprisingly well if we replace equation (|10|) with the 
modified formula 

M2 M|,(n) ^ M|,(ir-n-l) ^ 



K-1 n K -n-l 

In other words, some fermionic states at resolution K in the spectrum appear to be well 
approximated as two |-F1) particles with a combined momentum of K — 1. There is 
still one unit of momentum unaccounted for, and one proposal [|| is to assume that this 
momentum is associated with a topological configuration of the light-cone vacuum that 
makes it distinct from the trivial light-cone vacuum. In fact, in order to get a fermion 
state overall, we need to assume that this topological vacuum carries fermion quantum 
numbers. 

Although it is indeed possible to construct topologically distinct vacuum sectors in 
an SU(A^) light-cone gauge theory with the properties cited above |]13], [1^, the empirical 
success of equation ([TT| ) lends itself to an alternative interpretation. Firstly, we know 
the theory has a supersymmetric point for a very specific value of the fermion mass [^: 
TV? = g'^N/iT. As we reduce the fermion mass m to zero, one can show numerically that 
the mass degeneracy between the lightest boson-fermion supersymmetric pair is broken 
linearly as a function of th? / g'^N P|. However, if we consider a heavier supersymmetric 
pair, this behavior is not necessarily observed. For example, at the supersymmetric 
point {m? = q'^N/tt), the mass of the next-to-lightest fermion in the K = 15 DLCQ 
spectrum is approximately = 43.5, while its super-partner is identified in the K = lA 
DLCQ spectrum as a boson with massQ = 40.7 (in units g'^N/ir). The boson state 
is essentially a mixture of two-bit and four-bit single-trace states, while the fermion is 
observed to be essentially a mixture of three-bit single-trace states. 

If we carefully track these states as the fermion mass m is decreased to zero, we find 
in the massless limit that the fermion state has mass = 20.86, while its (former) 
^The boson and fermion masses are expected to converge together in the continuum hmit X — > oo. 
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superpartner has the new mass = 21.63. Surprisingly, the masses are still approxi- 
mately degenerate, so supersymmetry appears to be very weakly broken for these states. 
This suggests that the theory is in fact asymptotically supersymmetric (i.e. the spectrum 
of highly excited bound states is almost supersymmetric for any fermion mass m), an 
observation pointed out by Boorstein and Kutasov . One can use formula (|1^) to show 
that the boson is in fact a composite state of two freely interacting \F1) particles [Q. 
We conclude therefore that the mass of particular fermionic states at DLCQ resolution 
K may be estimated by an application of equation ( p!OD with K replaced with K — 1, 
since the composite boson at resolution — 1 is still approximately equal in mass with 
the fermion state at resolution K. This scheme provides a natural derivation of equation 
(0) from (|10D. 

At this point, we have made no reference to composite particles made from three 
or more single-particle states; to further validate our interpretation of the spectrum, we 
look for states of the form |F1) ® \F1) ® \F1). The mass Mp-^^f^^p-^^^p-^ is now given by a 
simple generalization of equation (p!OD: 



K - K, ^ K, ^ K, ' ^''> 
where Ki + K2 + K3 = K. For example, setting Ki = 3,K2 = 5, and K3 = 7 gives 
"^Fi(g)FiiS)Fi — 48.6915. This is indeed observed in the DLCQ fermionic spectrum at 
K = 15, and we observe an exact two- fold degeneracy at any finite A^. These states 
involve significant contributions from multi-trace Fock states, and are therefore only seen 
in a finite N analysis. 



4 Discussion 



In this work we have implemented DLCQ to study the spectrum of QCD^^^,]^ coupled 
to a massless Majorana fermion transforming in the adjoint representation of the gauge 
group SU(A^), where N is any integer > 3. We observed in Tables |1| and ^ that the 
masses of the stable "single-particle" states in the spectrum are very weakly dependent 
on A^. Evidently, considerable cancellations must be responsible for protecting masses 
from large 1/A^ corrections. It would be desirable to have a better understanding of this. 

After comparing predicted masses for the composite particle \F1)®\F1) (Table|^) with 
the actual DLCQ spectrum for bosons (Table ^), we observed that the finite A^ spectrum 
is consistently described as multi-particles with freely interacting stable constituents. 
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Therefore, the existence of a continuous spectrum above the threshold mass = AMp^ 
appears to be vahd not only for = oo but also for any finite N. This behavior 
is compatible with the viewpoint that the finite N theory is in a screening phase 0. 
Clearly, probing larger values of K will help clarify this issue, and we leave this challenge 
for future work. 

We also attempted to explain the fermionic multi-particle spectrum by assuming 
the theory was asymptotically supersymmetric; i.e. excited states in the spectrum are 
almost supersymmetric for any fermion mass, including m = 0. An explicit example 
was given. This enables one to connect the mass of a (sufficiently heavy) fermion with 
its corresponding bosonic 'super-partner'. Because the DLCQ resolution of a boson and 
fermion state differ by one unit of momentum (after imposing anti-periodic boundary 
conditions for fermions), we were able to derive naturally formula ( pAj) from equation 
(p!0|) describing the mass of two free particles. 

Finally, we comment that the finite N analysis we have adopted here affords inter- 
esting insights into large N theories in general; solving a theory for increasing values of 

enables one to study the large A^ limit in a way that is not accessible in the N = oo - 
or planar - approximation. In particular, subtleties can arise in the large A^ limit where 
single and multi-trace states co-exist on an equal footing. This is especially true for 
theories that may be susceptible to a phase transition - such as the model studied here 
- and may play a crucial role in understanding the fundamental degrees of freedom of 
large A^ gauge theories. 
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